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ABSTRACT 



Context. Gas-rich dwarf galaxies are probably the closest counterparts to primeval objects we can find in the local Universe, therefore 
it is interesting to study their evolution in different astrophysical contexts. 

Aims. We study the effects of interstellar clouds on the dynamical and chemical evolution of gas-rich dwarf galaxies. In particular, we 
focus on two model galaxies similar to IZwl8 and NGC1569 in comparison to models in which a smooth initial distribution of gas is 
assumed. 

Methods. We use a 2-D hydrodynamical code coupled with a series of routines able to trace the chemical products of SNell, SNela 
and intermediate-mass stars. Clouds are simulated by adding overdense regions in the computational grid, whose locations are chosen 
randomly and whose density profiles match observed ones. We consider both cloud complexes put at the beginning of the simulation 
and a mechanism for continuous cloud formation. The clouds are inherently dynamically coupled to the diffuse gas, and they experi- 
ence heat conduction from a hot surrounding gas. 

Results. Due to dynamical processes and thermal evaporation, the clouds survive only a few tens of Myr. Due to the additional cooling 
agent, the internal energy of cloudy models is typically reduced by 20 - 40% compared with models of diffuse gas alone. The clouds 
delay the development of large-scale outflows by mass loading, therefore helping to retain a larger amount of gas inside the galaxy. 
However, especially in models with continuous creation of infalling clouds, their bullet effect can pierce the expanding supershell 
and create holes through which the superbubble can vent freshly produced metals. Moreover, assuming a pristine chemical compo- 
sition for the clouds, their interaction with the superbubble dilutes the gas, reducing the metallicity. The resulting final metallicity is 
therefore generally lower (by ~ 0.2 - 0.4 dex) than the one attained by diffuse models. 
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1. Introduction 

Gas-rich dwarf galaxies are commonly classified into low 
surface-brightness dwarfs, called dwarf irregulars (dlrrs), and 
higher surface-brightness objects, usually called blue compact 
dwarf (BCD) galaxies. These classes of galaxies tend to have low 
metallicities, blue colors and complex and chaotic gas phases. A 
large fraction of these galaxies shows an ongoing star formation 
(SF) or at least hints that this process has been quenched in the 
recent past. In this case, these objects are commonly referred to 
as starburst galaxies and their gas consumption timescales are 
much shorter than the Hubble time (Kennicutt 1998), making 
this a transient phase of their evolution. 

Owing to the energy released by stellar winds and super- 
novae (SNe), intense episodes of SF are also associated to the 
development of galactic winds or at least of large-scale outflows. 
The broad distinction between these two phenomena is the final 
fate of the outwards-directed flow of gas: galactic winds gen- 
erally exceed the escape velocity while outflows do not, there- 
fore they tend to recede towards the center of the galaxy. Clear 
signatures of outflows are present in NGC1705 (Hensler et al. 
1998, He ckman et al. |200TT ), NGC1569 (Marti n, Kobulnicky & 
Heckman|2002]), NGC3079 (Cecil et al. 120011 IZwl8 (Martin 
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1996), NGC3628 (Irwin & Sofue 1996) among others. Perhaps 
the best examples of large-scale outflows driven by SN feedback 
are at large redshifts (Pettini et al. 119981 Pettini et al. I2001I I. 
Although it is not certain, in any of the above-mentioned ob- 
jects, that the metals will definitely leave the parent galaxy, in- 
direct hints of the ubiquity of galactic winds are given by the 
mass-metallicity relation (Tremonti et al. 2004; Davee, Finlator 
& Oppenheimer 2006 ) and effective yields (Garnett 2002). 

From a theoretical point of view, the study of the evolu- 
tion of gas-rich dwarf galaxies through numerical simulations 
has been performed by several authors in the recent past. The 
overall picture is that the occurrence of large-scale outflows 
is initially driven by the thermal pressure of a very hot, high 
pressurized gas and is favored by a flat distribution of the in- 
terstellar medium (ISM), which allows an easy vertical trans- 
port of material. However, since the transport of gas along the 
disk is very limited, outflows are not able to eject a significant 
fraction of the ISM, whereas the fraction of ejected metals can 
be very large (D'Ercole & Brighenti 1999; MacLow & Ferrara 
[T9991 Recchi, Matteucci & D'E rcole l^OOTl hereafter RMD). For 
NGC1569, Martin et al. (2002) derived a supersolar metal con- 
tent in the galactic wind from X-ray spectra but also advocated 
mass-loading of it with the ISM. 

Most of these studies, however, have focused on flows in 
homogeneous media, neglecting the multiphase nature of the 
ISM, although several attempts to perform multiphase hydro- 
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dynamical simulations have been made in the past, particularly 
using the so-called chemodynamical approach (Theis, Burkert 
& Hensle r [19921 Rieschick & Hensler |2000l Hensler, Theis & 
Gallagher 2004). The multiphase nature of the ISM, in partic- 
ular its dumpiness, is observationally well established in dwarf 
galaxies (Cecil et al. l200"Tl Cannon et al. 120051 Leroy et al. [2006b 
and it has a solid theoretical background with the seminal work 
of McKee & Ostriker dl977l l. According to this model, the ISM 
is composed by a cold neutral phase (representing the cores of 
molecular clouds), confined by a warm medium (with temper- 
atures of the order of 10 4 K) and these two phases (which are 
in pressure equilibrium) are embedded in a hot, diluted inter- 
cloud medium (HIM), continuously produced by SN explosions 
and stellar winds. Sufficiently dense clouds can pierce the HIM 
without being swept up, so they can become embedded therein 
(Vieser & Hensler |2007bl ). At the interface between clouds and 
HIM, condensation-evaporation processes establish the final fate 
of the cloud and its impact on the development of a galactic 
wind. 

In two previous papers, we have studied the dynamical and 
chemical evolution of model galaxies similar to IZwl8 (Recchi 
et al. 120041 hereafter Paper I) and NGC1569 (Recchi et al. l2006l 
hereafter Paper II). The main results can be briefly summarized 
as follows: 

- most of the analyzed models develop large-scale outflows. 
These outflows carry out of the galaxy mostly the chemical 
elements freshly produced during the most recent episodes 
of SF, with large escape fraction of metals with delayed pro- 
duction (like Fe and N). 

- Models with very short burst(s) of SF can cool and mix 
the newly formed metals in a very short timescale, whereas, 
when the SF is more complex, most of the metals are either 
directly ejected outside the galaxy through galactic winds or 
are confined in a too hot medium, therefore cannot contribute 
to the chemical enrichment of the warm ionized medium ob- 
served by emission lines from the H n gas. 

- Models with complex and long-lasting SF episodes repro- 
duce the chemical composition and the abundance ratios of 
the above-mentioned galaxies much better than models with 
bursting SF. 

In this paper we simulate models with structural parameters 
similar to IZwl8 and NGC1569. We increase arbitrarily the gas 
density of some specific regions of the computational grid, in 
order to create a "cloudy" phase, and we address the question 
how and to which extent a "cloudy gas phase" alters the former 
results. The clouds possess a specific density profile and can be 
either added at the beginning of the simulation or continuously 
created during the evolution of the model. We then analyze the 
differences between the dynamical and chemical evolution of 
these models with the ones presented in Paper I and Paper II. We 
point out that, at variance with the above-mentioned works, in 
this paper we will not specifically look for the best initial setups 
and the best assumptions in order to reproduce chemical and dy- 
namical features of well-known objects. We will just stress the 
main variations produced by a clumpy initial setup. For this rea- 
son, we will also consider models which failed in Paper I and II 
at reproducing the observations of IZwl8 and NGC1569. 

The paper is organized as follows: in Sect.|2]we briefly recall 
the evolution of a cloud embedded in a hot medium; in Sect. [3] 
we present the model and the adopted assumptions in the simula- 
tions. Results are presented in Sect. |4](models with clouds fixed 
at the beginning of the simulation) and in Sect. [5] (continuous 



creation of clouds). Finally, a discussion and some conclusions 
are drawn in Sect. [6] 

2. The dynamics of clouds embedded in a hot 
phase 

The ubiquitous coexistence of cool to warm clouds in the hot 
phase of the ISM has attracted during the recent years broad 
attention on the interaction of such clouds with this tenuous 
and hot medium. The studies deal with three major effects: the 
influence of heat conduction, of shock fronts, and of the dy- 
namics of the flowing hot gas, like e.g. with or without the 
presence of large-scale galactic outflows, on the evolution of 
clouds (Ha rtquist et al. 119861 Murray et al. 119931 Ferrara & 
Shche kinov |19931 Vi etri, F errara & Miniati IT9971 F ragile et al. 
120041 Marcolini et al. |20031 Vieser & Hensler |2007bl among oth- 
ers). Despite the very large variety of adopted methodologies, 
astrophysical contexts and involved physical processes, and in 
spite of clearly defined problems, a broad variety but not yet 
uniqueness of issues can be summarized as follows: 

- Moving from the idealized situation treated in analytical 
thermal conduction studies (e.g. Cowie & McKee 119771 ) to 
saturated heat conduction and self-gravitating clouds can 
change the results from evaporation to condensation for the 
same state of the hot gas and the same cloud mass model 
(Vieser & Hensler l2007a1 ). 

- A clumpy medium embedded in a hot flow produces mass 
loading, namely the seeding of material, ablated from the 
clouds, into the global flow. It has been demonstrated that 
this kind of phenomena helps in clarifying the X-ray emis- 
sion in starburst galaxies like M82 (Suchkov et al. 119961 
Strickland & Stevens 120001 Marcolini et al. l2005l 

- A single cloud overrun by a shock wave can be crushed 
within the so-called crushing time (i.e. the time needed for 
the internal forward shock to cross the cloud and reach its 
downstream surface) and will be destroyed to smaller frag- 
ments if cooling dominates. But vice versa it also evaporates 
in a parameter regime with exceeding thermal conduction 
(Orlando etal. EES). 

- In a complex of clouds, if the cloud separation transverse 
to the flow is smaller than some critical value (a few times 
the typical cloud radius, the exact value depending on the 
authors), the clouds will merge into a single structure before 
the hot flow destroys them. 

- Thermal conduction helps in stabilizing the surface of the 
cloud, making it less susceptible to Kelvin-Helmholtz and 
Rayleigh-Taylor instabilities (Orlando et al. 2005 1 Vieser & 
Hensler :2007b). It can also generate an inward-propagating 
shock wave able to compress the core of the cloud. 

In the present work we do not intend to simulate in great de- 
tail the interaction of a cloud or of a cloud complex with a diffuse 
hot medium, as made by the previously cited authors. We in- 
stead simulate galaxy models similar to well observed and stud- 
ied gas-rich dwarf galaxies, relaxing the hypothesis of a smooth 
initial gaseous distribution (as assumed in Paper I and Paper II) 
and analyze how the inclusion of a clumpy medium changes the 
thermal and chemical evolution of the ISM. Indeed, the reso- 
lution required to properly take into consideration conductive 
fronts surrounding clouds is of the order of 0. 1 pc (Marcolini et 
al. 120051 Vieser & Hensler 2007a,b), extremely computationally 
demanding in a simulation in which the large-scale evolution of 
the galaxy, up to distances of several kpc has to be taken into 
account. 
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3. Model description 

3. 1 . The numerical code 

The simulations are performed by means of a 2-D hydrocode in 
cylindrical coordinates based on a second-order upwind scheme 
(Bedogni & D'Ercole 1986). The hydro solver is coupled with 
routines able to follow in detail the chemical and dynamical 
feedback on the galaxy as a consequence of SNell, SNela and 
winds from intermediate-mass stars (IMS). 

The chemical evolution has also an impact on the dynamics 
of the system, due to the assumption of a metallicity-dependent 
cooling function (Bohringer & Henslerfl9891). We point out that, 
as in the previous papers, when we plot abundances or abun- 
dance ratios produced by our models, we exclude grid points at 
temperature above 2 ■ 10 4 K. This is because the gas at these 
temperatures would be undetectable with optical spectroscopy 
and its metallicity could be guessed only through X-ray analysis, 
whose use in dwarf galaxies is still quite uncertain (Martin et al. 
120021 Ott, Walter & Brinks 120051 The code has been described 
in detail in RMD and newer implementations and improvements 
are reported in Paper I and Paper II, therefore we refer the read- 
ers to these papers for technicalities. 

Given the importance of thermal conduction on the shap- 
ing and final fate of clouds embedded in a hot medium, we just 
briefly recall the numerical method adopted to treat this physical 
phenomenon in our code. To solve the heat transport equation, 
the operator splitting method is adopted and the one-dimensional 
problem is solved through the Crank-Nicholson method (see also 
D'Ercole & Brighenti 1999). A saturated heat flux (Cowie & 
McKee 119771 ) is adopted if the mean free path of electrons is 
larger than the temperature scaleheight. 

3.2. Model parameters 

The initial configurations of our models are aimed at repro- 
ducing the main structural parameters of two very well stud- 
ied gas-rich dwarf galaxies: IZwl8 and NGC1569. The initial 
setup is taken from our previous models and described in RMD 
(for IZwl8) and in Paper II (for NGC1569). As described in 
Sect. 13.31 the gaseous distribution is made clumpy, either per- 
turbing the initial setup or adding clouds as the galaxy evolves, 
at the same rate as the SF. We also vary the initial mass function 
(IMF) of the stars. The models are identified through the nota- 
tion XYZW, where X refers to the initial setup (I: setup similar 
to IZwl8; N: setup similar to NGC1569), Y takes into consider- 
ation whether the clouds are put in the initial setup of the galaxy 
or are continuously created (B: clouds present from the begin- 
ning; C: continuous creation of clouds). The third index refers to 
the adopted IMF: S is for the Salpeter index (x=1.35). A flatter- 
than-Salpeter index is also tested: A (x=0.95) and a steeper index 
(x=1.7) is denoted with K. 

Finally, also the yields from IMS are allowed to change (R 
for Renzini & Voli 1981, V for van den Hoek & Groenewegen 
119971 and M for Meynet & Maeder 120021 ). For instance, the 
model called IBSR starts with a setup aimed at simulating 
IZwl8, puts the clouds at the beginning of the simulation, as- 
sumes the yields of Renzini & Voli ( 198 1) for IMS and a Salpeter 
IMF. The yields from massive stars are taken from Woosley & 
Weaver ( 1995), unlike for model M, where, for self-consistency 
reasons, the yields of Meynet & Maeder (2002} are assumed. 
In this case, the upper mass is also 60 M , at variance with 
the 40 M value adopted for Woosley & Weaver yields. For the 
SNela, the formulation of Matteucci & Recchi (2001) has been 



adopted, with nucleosynthetic yields taken from the model W7 
of Nomoto, Thielemann & Yokoi ( 1984). The model parameters 
are summarized in Table Q] 

As applied in Paper I and Paper II, the SF is gasping (i.e. long 
episodes of SF are separated by short periods of quiescence). 
More details about the SF history of various models are given in 
the corresponding sections. We will be referring to diffuse mod- 
els any time we will consider a model, similar to the described 
cloudy one, but with a homogeneous initial density distribution, 
described in detail either in Paper I or in Paper II. 

3.3. Cloud description 

A random generator identifies grid points in the galactic region 
which are the cores of the clouds. In the models in which clouds 
are put "ab initio" in the setup of the model (models identified 
with the second index "B", see Sect. 13.21 ). the central density is 
decided a priori (100 cm -3 ), and the total number of clouds is 
25. The exact number of clouds does not play a significant role; 
tests have been performed also with 50 or 75 clouds (reducing 
their masses) and this does not affect significantly the results. 
The probability of finding a cloud in a particular grid point is 
proportional to the gas density. Due to the flattened distribution 
of the ISM at the beginning of the simulation (see RMD or Paper 
II), this gives also a larger probability for the clouds to be found 
close to the disk. The density of the grid points outside the clouds 
is then reduced in order to get a consistent final gaseous mass. 
In models with continuous creation of clouds (identified with the 
second index "C") the mass is decided a priori: it is assumed that 
the clouds are created at the same rate as the SF rate. 

As explained in Paper I, all the stars formed within an inter- 
val of time At are treated as a single stellar population. After this 
interval of time (typically 10 6 yr), a cloud is created, whose mass 
equals the mass of the single stellar population. Moreover, these 
clouds are given an initial infall velocity of 10 km s _1 . For any of 
these models, the clouds are approximately shaped according to 
a radial density profile p c \ oc R~ li7 , where R c \ is the distance from 
the center of the cloud and the exponent of the power-law den- 
sity profile is taken from observations (de Heij, Braun & Burton 
120021 Churchill, Vogt & Charlton|2002 Tatematsu et al. [2004l) . 
It is important to stress that, given the dimensionality and the 
symmetry of our numerical code, the clouds are not spherical, 
but ring-like structures. The clouds are assumed to have primor- 
dial chemical composition (i.e. no metals in them) and the tem- 
perature is set to 10 3 K. Since this is also the minimum temper- 
ature allowed for the gas (see Paper I and Paper II), the clouds 
are not in pressure equilibrium with the surrounding medium. In 
Sect l4. 1 . il we offer arguments demonstrating that this assump- 
tion is not unrealistic and we briefly describe models in which 
this assumption is relaxed and clouds are put in pressure equilib- 
rium with the surrounding medium. 

4. Results of models with fixed clouds. 

4.1. Model IBSR 

4.1.1. Dynamical evolution 

We use as a prototypical model with a fixed complex of clouds a 
setup similar to the model SR2 analyzed in Paper I. It has there- 
fore the same SF history adopted therein, namely a long, moder- 
ate episode of SF lasting 270 Myr, a quiescent period of 10 Myr 
and a more vigorous burst (5 times more intense than the first, 
long-lasting SF episode) lasting 5 Myr (Aloisi, Tosi & Greggio 
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Table 1. Summary of model parameters 



1V1UQC1 


Setup 


cloud 
creation 


TIx/Th c 1 nnp v 
llvir SlUJJC A 


TMS viplrk* 


IBSR 


IZwl8 


no 


1.35 


RV81 


IBSV 


IZwl8 


no 


1.35 


VG97 


IBAV 


IZwl8 


no 


0.95 


VG97 


ICSV 


IZwl8 


yes 


1.35 


VG97 


ICKV 


IZwl8 


yes 


1.70 


VG97 


NCSM 


NGC1569 


yes 


1.35 


MM02 



RV81: Renzini & Voli ( HMB ; VG97: van den Hoek & Groenewegen ( fl997l >; MM02: Meynet & Maeder (2002 1. 



1999 ). We add randomly 25 clouds, according to the procedure 
described in Sect. 13.31 The initial gas distribution is shown in 
Fig. [T] It is worth reminding that, since the central density of 
the clouds is decided a priori, each cloud has a different mass, 
the lightest ones being close to the center. In Fig. [2] we plot the 
logarithmic distribution of the cloud masses. As we can see, the 
clouds span almost 2 orders of magnitude in mass, ranging from 
~ 10 3 5 to ~ 10 5 5 M , with a mean value of 6.26 • 10 4 M . 



Density 




R (Kpc) 

Fig. 1. Density contours for model IBSR at the beginning of the 
simulation. The density scale (in g cm 4 ) is on the right-hand 
strip. 



The evolution of this model in the first ~ 120 Myr is shown 
in Fig. [3] This model is able to develop a large-scale outflow in 
a timescale of the order of ~ 75 Myr, only slightly delayed com- 
pared to the diffuse model SR2 (for comparison, see figs. 1 and 
2 of Paper I). This model shows however a much more distorted 
density structure and much larger eddies and regions of ther- 
mal instabilities. As we will see better later on, the clouds are 
destroyed in a relatively short timescale, but nevertheless they 
leave an imprint on the development of the outflow, strongly in- 
fluenced by the patchiness of the medium in the central region 
of the galaxy. 

To quantify this effect, we calculate the total thermal energy 
of model IBSR during the first ~ 100 Myr and compare it with 
the value found for the model SR2 in Paper I. The thermal en- 
ergy is calculated inside a region R < 1 kpc and < 730 pc, which 
we have called 'galactic region' in RMD and which nearly coin- 
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Fig. 2. Histogram of the logarithmic distribution of the cloud 
masses (in M ) for the model IBSR. The dashed line represents 
the mean value (6.26 ■ 10 4 M ). 



cides with the region where stars are distributed. This compari- 
son is shown in Fig. [4] As we can see, the total thermal energy 
budget is clearly affected by the presence of the clouds, leading 
to a reduction of -20% of the thermal energy deposited into the 
system, i.e. it leads to a ~ 20% increase of the radiative losses. 
Incidentally we notice that this value of the total thermal en- 
ergy (corresponding to the explosion energy of just a few SNe) 
is strongly affected, more than from the radiative cooling of the 
superbubble, from our assumption of a low thermalization ef- 
ficiency of SNell. An extensive discussion about this debated 
parameter can be found in RMD and Paper I. 

In Fig.|5]we also show the average radius of the superbubble 
for the two above-mentioned models. The cloudy model IBSR 
allows at the beginning a slightly faster expansion of the super- 
bubble. This is due to the fact that, in order to reproduce the same 
total mass, in a cloudy model the density of the diffuse medium 
has to be reduced. Moreover, the presence of clouds strongly dis- 
torts the shape of the supershell and the highly pressurized gas 
inside the cavity can more easily find regions of lower pressure, 
from which it is easy to pierce the shell and break out. This cre- 
ates the tongues visible in Fig. [3] As time goes by, the larger 
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0.5 1 0.5 1 0.5 1 0.5 1 
R (Kpc) R (Kpc) R (Kpc) R (Kpc) 

Fig. 3. Density contours and velocity fields for model IBSR at four different epochs (evolutionary times are labeled in the box on 
top of each panel). The logarithmic density scale (in g cirT 3 ) is given in the strip on top of the figure. In order to avoid confusion, 
velocities with values lower than 1/10 of the maximum value (indicated for each panel in the upper right box) are not drawn. This 
is valid also for Fig. [6] 



radiative losses experienced by the model IBSR slow down the 
expansion of the superbubble and at ~ 50 Myr the average su- 
perbubble radius of the diffuse model overcomes the one of the 
cloudy medium. 

We can better analyze the influence of clouds on the early 
development of galactic winds by zooming in the central re- 
gion of the computational grid during the first tens of Myr. 
A plot of the density profiles, velocity fields and temperature 
profiles of the model IBSR in the first ~ 20 Myr is shown in 
Fig. [6] We can first notice in this plot that the clouds, even 
if they are not encompassed in the superbubble, tend to ex- 
pand. This is because they are not put in pressure equilibrium 
with the surrounding ISM, therefore their lifetime is relatively 
short (a few tens of Myr), irrespective of the presence of a sur- 
rounding HIM. We have also run models in which the clouds 
are put in pressure equilibrium with the surrounding ISM. This 
has been obtained by simply reducing the temperature of the 
clouds up to the value which compensate for the pressure of 



the ISM at the border of the cloud. This model shows the same 
behavior of the other models once the clouds are encompassed 
within the superbubble cavity, namely they are quickly evapo- 
rated. The clouds outside the superbubble show of course a dif- 
ferent behavior, being more stable than the previously described 
ones, but this does not affect the global evolution of the model. 
Moreover, growing evidences are accumulating, both observa- 
tionally (Ballesteros-Paredes, Hartmann & Vazquez-Semadeni 
119991 Hartmann, Ballesteros-Paredes & Bergin [200U Hartmann 
120031 and theoretically (Elmgreen [20001 Heitsch et al. 120061) 
that molecular clouds are transient structures rather than well 
defined objects in quasi-equilibrium states. Therefore they must 
have a relatively short lifetime, as in our simulated clouds. Given 
the similarities between the behavior of equilibrium and non- 
equilibrium clouds and due to the more complex and slow com- 
putation of equilibrium clouds, we focus from now on only on 
models in which the clouds are not in pressure equilibrium. 
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Fig. 4. Thermal energy (in units of 10 erg) for the model IBSR 
(solid line) and for the reference diffuse model SR2 (Paper I) 
(dashed line). 



Fig. 5. Average superbubble radius (in kpc) for the model IBSR 
(solid line) and for the reference diffuse model SR2 (dashed 
line). 



The clouds create a very patchy temperature distribution in 
the first 10-12 Myr but the evaporation process continues up 
to the moment (after 20 Myr) where the temperature inside the 
superbubble is almost uniform. We can notice once more a su- 
pershell which, owing to the interaction with clouds, strongly 
deviates from the spherical geometry, allowing therefore larger 
shears and eddies (and consequently a reduced thermal energy 
of the bubble), but also tongues and fingers from which it is eas- 
ier to leak out of the superbubble. This is the reason why, after 
~ 100 Myr some outflowing gas is already ~ 3 kpc above the 
galactic disk (third panel of Fig. |3]l although the average super- 
bubble radius is at this stage smaller than for the diffuse model 
SR2(Fig.|5}. 

4.1.2. Chemical evolution 

From a chemical point of view, the expected effect of clouds 
is to dilute the hot metal-rich gas through evaporation of the 
(metal-poor) clouds, allowing for a reduction of the metallic- 
ity without altering the abundance ratios (Koppen & Hensler 
2005). However, the reduced thermal energy and, consequently, 
the reduced escape fraction of metals from the galactic region 
(Sect. I4.1.1I ) should produce an increase of the metallicity of the 
galactic ISM. In this case, the abundance ratios are affected if 
the ejection efficiencies depend on the different chemical species 
(here we simply define ejection efficiency as the fraction of met- 
als outside the galactic region compared to the total amount 
which has been synthesized). To disentangle these two compet- 
ing effects, we analyze the differences in the chemical evolution 
between diffuse and cloudy models. 

The comparison of 12+log(0/H), log(QO), log(N/0) and 
[O/Fe] of the warm ionized phase is presented in Fig. [7] In this 
plot we also show the evolution of a prototypical model of con- 



tinuous creation of cloud (model ICSV, see Sect. I5.1.2t . The 
evolution of a-elements is not significantly altered by the pres- 
ence of clouds (the difference being always around ~ 0. 1 dex), 
but we can notice less negligible differences (of the order of 0.2 
dex) in the log(N/0) abundance ratio. Indeed the diffuse model, 
due to the reduced radiative losses, attains a larger fraction of 
metals with temperature above 2 • 10 4 K, therefore excluded by 
this plot. Nitrogen is mostly produced during the thermal puls- 
ing phase by AGB stars of masses ~ 4 - 7 M therefore with a 
delay compared to the prompt production of oxygen. Soon after 
its production, nitrogen is also located inside a hot cavity (carved 
by SNe) most likely than oxygen. Moreover, as demonstrated in 
RMD, the ejection efficiency of nitrogen can be larger than the 
one of a-elements, again favoring the decrease of N/O in the dif- 
fuse model, for which the development of a large-scale outflow 
is anticipated. 

Incidentally, we can notice that the assumption of a cloudy 
medium worsens the agreement between the predicted log (N/O) 
and the observations. However, as we have pointed out in Paper 
I, the Renzini & Voli dl98H yields tend to overestimate the ni- 
trogen production and only the assumption of the Meynet & 
Maeder (2002) set of yields can reconcile the prediction of the 
models with the observed abundance ratios in IZwl8. We stress 
once again that the main goal of this paper is a study of the effect 
of a cloud complex on the dynamical and chemical evolution of 
a gas-rich dwarf galaxy rather than the attempt to exactly repro- 
duce the chemical features of specific objects. 
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Fig. 6. Density contours and velocity fields (left panels) and 
temperature contours (right panels) for the central region of 
model IBSR at four different epochs (evolutionary times are la- 
beled in the box on top of each right panel). Logarithmic scales 
are given on top of each column of panels. 



4.2. Model IBSV 
4.2.1. Dynamical results 

The model IBSV differs from the model IBSR only on the as- 
sumed IMS yields (van den Hoek & Groenewegen |1997"l instead 
of Renzini & Voli 1981; see Table[T|l. Dynamically, the only dif- 
ference is therefore a variation of the cooling rates, due to our 
assumed metallicity-dependent cooling function. It is interesting 
to quantify this effect through a direct comparison of models dif- 
fering solely in the adopted nucleosynthetic yields. This compar- 
ison is shown in Fig.[8]at two evolutionary times: 50 Myr (lower 
panels) and 100 Myr (upper panels). Overall, the agreement be- 
tween the dynamics of these two models is very good. However, 
the extremely non-linear character of superbubble evolution can 
be noticed in this plot and small differences in the energy budget 
of the model (due to different cooling rates) can produce signifi- 
cant changes in the dynamics. In particular, in model IBSR after 
100 Myr some gas is already flowing out of the galaxy through 
a narrow nozzle, whereas in model IBSV the occurrence of a 
large-scale outflow is slightly delayed. At this point, the total 
mass of gas within the galaxy is ~ 9% larger in model IBSV, a 
small but non-negligible factor. 




100 200 300 

I (Myr) 



Fig. 7. Evolution of 12 + log (O/H) (top panel), log (C/O) (sec- 
ond panel), log (N/O) (third panel) and [O/Fe] (bottom panel) for 
a prototypical model with a fixed initial cloud complex (IBSR 
model) (solid line) and a prototypical model with continuous 
creation of cloud (ICSV model) (dashed line). These models are 
compared with a model of similar mass but smooth ISM distri- 
bution (dotted line). The superimposed shaded areas indicate the 
observed values found in literature (if available), with relative 
error-bars. 



4.2.2. Chemical results 

The comparison of log(C/0) and log(N/0) in models IBSR and 
IBSV is shown in Fig. [9] In this plot we do not show the evo- 
lution of oxygen because its production in IMS is negligible. 
Since the prescriptions for the yields from massive stars are 
the same (Woosley & Weaver [19951 1, we do not see significant 
differences in the two models. The yields of van den Hoek & 
Groenewegen (1997) produce substantially less carbon and ni- 
trogen compared to Renzini & Voli (1981) yields (Chiappini, 
Romano & Matteucci 2003). The difference is particularly sig- 
nificant for what concerns N, whose production is approximately 
halved, resulting therefore in a N/O ~ 0.3 dex lower than in the 
model IBSR. The same difference has been produced by diffuse 
models with different IMS yields (e.g. Recchi et al. 120021 Paper 
I), indicating that the different dynamics of the cloudy model do 
not significantly affect this behavior. 

4.3. Varying IMF: model IBAV 

As shown in Table Q] a model similar to IBSV but with flatter 
IMF is also considered. This models has the same SF history 
considered so far (i.e. the one derived from the work of Aloisi 
et al. 1 19991 1, but the energy injection rate is much larger. To 
be more precise, the total energy release is ~ 2.5 times larger 
than in IBSV. Since the gas binding energy remains the same, in 
spite of the larger radiative losses due to the interactions clouds- 
HIM, this energy release is enough to unbind all the gas initially 
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Fig. 8. Density contours for model IBSR (left panels) and model IBSV (right panels) at two evolutionary times: 50 Myr (lower 
panels) and 100 Myr (upper panels). The density scale (in g cm 4 ) is on the right-hand strip. 



present in the galactic region (an ellipsoid of dimensions ~ 1000 
x 730 pc, see RMD) in ~ 250 Myr. This complete blow-away 
does not happen in IBSV, where a large-scale outflow occurs in 
the polar direction but most of the gas close to the disk remains 
in the galactic region. 

5. Results of models with continuously created 
clouds. 



As described in Sect. 13.31 in this set of models we produce a 
cloud each At yr (typical value 10 6 yr), having a mass equal to 
the total amount of gas turned into stars in the same interval of 
time. In the framework of simple closed-box models of chemi- 
cal evolution, this case, (infall rate equal to the SF rate) is called 
extreme infall (Larson [1972| l and leads to the simple expression 



for the metallicity Z = yz[l — e 



], where // is the gas mass 



fraction and yz is the total yield (i.e. the ratio between the total 
mass in metals newly formed and the amount of mass locked up 
in low mass stars and remnants). The relaxation of the instanta- 
neous recycling approximation and the inclusion of dynamical 
effects (winds and mixing and cooling of metals) changes this 
finding but Koppen & Edmunds (1999) demonstrated that the 
ratio between infall and SF rate is the determining factor in the 
chemical evolution of galaxies. The clouds are given an infall ve- 
locity of 10 km s along the polar direction, their location in the 
computational grid is randomly chosen and their profile is again 
Pel x ^Ci 1 ' 7 ' ^ ut ^ n this case the central density is constrained by 
the total mass of the cloud and by its location rather than being 
constant for each cloud. 



5.1. Model ICSV 
5.1.1. Dynamical results 

As prototype of this group of models, we use a setup similar 
to IBSV, the only difference being the mechanism of cloud for- 
mation. Given the assumed SF history, during the first episode 
the clouds have a mass of 6 • 10 3 M , which increases to 3 
• 10 4 M during the last burst. Snapshots of the evolution of 
this model in the first ~ 55 Myr are shown in Fig. [TO] In this 
figure the shocks created by the clouds in their descent to the 
galactic disk are quite evident, in particular in the bottom row 
of panels. In particular, a bow shock is created around the cloud 
and a reverse shock is generated downstream behind the cloud, 
leaving an underdensity region behind it. The structure is also 
highly Kelvin-Helmholtz unstable. The timescale for the growth 
of Kelvin-Helmholtz instabilities is approximately 



0.5 



(1) 



inf 



(Chandrasekhar |T961b where q is the ratio between the cloud 
and the intercloud densities, v; n f is the infall velocity of the cloud 
(relative to the local ISM) and k is the wavenumber of the unsta- 
ble mode. The most unstable models are the ones with k ~ r~ l 
(where r c is the radius of the cloud), leading to a ?k-h between 
10 and 20 Myr with our parameters (depending on the size of 
the clouds, which is not constant). This is therefore also the 
timescale for the fragmentation of the cloud and its mixing with 
the local ISM. At later times, a non-negligible probability exists, 
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Fig. 9. Evolution of log (C/O) (upper panel) and log (N/O) 
(lower panel) for models IBSR (solid line) and IBSV (dotted 
line). 



given our simplified assumptions, that a cloud is created directly 
inside the expanding superbubble. An example is visible in the 
upper right panel of Fig.[10](at (R, z) ~ (20, 200) pc). In this case, 
given the much larger density ratio q (between 10 4 and 10 5 ), the 
Kelvin-Helmholtz timescale becomes larger than the time con- 
sidered in our simulations (a few hundreds of Myr). The cloud 
is therefore ablated by the flow of gas pushed by the exploding 
SNe (creating mass loading) and evaporated by the large tem- 
perature of the cavity (few 10 6 K) in a timescale of the order of 
few tens of Myr. Occasionally, clouds are created close enough 
in space and time, such that mutual interaction between clouds 
manifests. In this case, the clouds form coherent structures like 
the ones described by Poludnenko, Frank & Blackman (2002) 
before being evaporated. 

Due to the combined effect of thermal evaporation and bow 
shocks, the thermal energy budget is still smaller (by ~ 15%) 
than for model IBSR (plotted in Fig. |4j. Owing to the larger 
radiative losses and to the ram pressure of the infalling clouds, 
a galactic wind develops after ~ 130 Myr, delayed compared to 
model IBSR. However, as pointed out in Sect. 14.1.11 the infall 
of clouds can help in structuring and fingering the supershell, 
creating ways out for the hot gas. Indeed, the ram pressure of 
the infalling clouds is p mm = pism ■ vjL ~ 10" 11 erg cm~ 3 , of 
the order of the ram pressure of the expanding supershell and 
larger than the thermal pressure of the hot cavity. Moreover, it is 
worth pointing out that, for this set of models, the total gaseous 
mass is kept constant by the assumption that the cloud creation 
rate should balance the SF rate (see Sect. 13 . 3b , at variance with 
model IBSR, in which a fraction of gas (at a rate of 6 • 10~ 3 M 
yr _1 , see Paper I) is continuously turned into stars. 



In Paper II (sect. 4.7) we have analyzed the effect of the in- 
fall (along the polar direction) of a very large and very massive 
cloud. In this case, the development of an outflow is completely 
hampered and all the metals freshly produced by the ongoing SF 
are trapped inside the galactic region, preventing the develop- 
ment of differential outflows, a natural outcome of this kind of 
simulations. 

As we have seen, in model ICSV the "cap" effect is less sig- 
nificant: it helps reducing the total thermal energy able to drive a 
large-scale outflow, but the clouds are dissolved on a timescale of 
the order of a few tens of Myr. Moreover, as seen in Sect. 15.1.11 
they can pierce the supershell, creating funnels for the free flow 
of the hot, high-pressure gas. They can therefore (slightly) de- 
lay the development of an outflow, but they cannot prevent it. 
Consequently, the chemical evolution of this model is strongly 
affected by the dilution effect of the clouds. In particular, the pro- 
cess of cloud dissolution described in Sect. 15.1.11 continuously 
allows a mixing of the metals with pristine gas. This effect re- 
sults to be, in our simulation, much larger than the "cap" effect of 
the infalling clouds. We can see the chemical evolution of model 
ICSV in Fig. [7] The final oxygen abundance is ~ 0.3 dex smaller 
that the oxygen attained by model IBSR. Even more important 
is the fact that, after ~ 90 Myr, the oxygen abundance mildly but 
constantly decreases as function of time. This is due to the fact 
that the selective loss of metals has not been suppressed and that 
the continuous creation (and subsequent disruption) of clouds 
mixes the ISM with unpolluted gas. Very significant is also the 
effect on the N/O abundance ratio (a difference of ~ 0.6 - 0.7 
dex), but we stress once again that in this case the main reason 
of this difference is the choice of IMS yields. As we can see from 
Fig.|9l the final log(N/0) of model IBSV is ~ 0.3 - 0.4 dex larger 
than model ICSV. 

We also tested a model with a larger initial total mass of gas 
at the beginning of the simulation. The setup of this model is 
similar of model SV3 of Paper I (i.e. a total initial mass of 3 x 
10 7 M instead of the standard value of 1.7 x 10 7 M Q ), but in 
it we apply the same procedure of continuous creation of clouds 
analyzed for model ICSV. The slightly lower density contrast 
parameter q does not significantly affect the overall cloud-ISM 
interaction process and the dissolution timescale of the clumps 
but, due to the larger ISM pressure, the development of a large- 
scale outflow is largely delayed. It occurs only ~ 250 Myr after 
the beginning of the SF process. 

5.2. Model ICKV 
5.2.1. Dynamical results 

This model has the same setup and same SF history of 
model ICSV, the only difference being a steeper (x=1.7) IMF. 
Consequently, the energy return rate is much smaller than the 
above-considered model (about a factor ~ 2) and, since the bind- 
ing energy of the gas is the same, the reduced power of the burst 
has deep consequences on the development of a large-scale out- 
flow. In this model, a break-out of the superbubble with conse- 
quent outflow of gas happens at around ~ 180 Myr, but its inten- 
sity is very mild and the continuous infall of clouds is sufficient 
to suppress its further development and to close the funnel. The 
final structure (after ~ 300 Myr) is an elongated ellipsoid of ~ 
700 x 220 pc, with just some traces of gas which has managed to 
leak out and flows freely, mainly along the polar direction. Later 
on, no SF is occurring anymore, therefore the process of cloud 
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Fig. 10. Density contours for the warm gas for model ICSV at 9 evolutionary times (labeled in Myr at the top right corner of each 
panel). The density scale (in g cm -3 ) is on the right-hand strip. 



formation is suppressed. Although the input of energy is con- 
tinuous in these models (due to the contribution of SNela), the 
produced energy is not enough to break-out again, therefore the 
supershell tends to recede towards the center of the galaxy (see 
Recchi & Hensler 2006). Only ~ 3% of the gas produced during 
the two episodes of SF leave the galactic region at the end of the 
simulation. 



5.2.2. Chemical results 

The chemical evolution is characterized by a strong bias to- 
wards low- and intermediate-mass stars, therefore its produc- 
tion of a-elements is strongly reduced. The oxygen abundance 
smoothly increases during the first ~ 180 Myr (approximately 
up to the break-out), then it stays almost constant. As stressed in 
Sect. 15.1.21 this is mainly due to the dilution effect of the clouds 
(as we have seen, the differential loss of metals is very limited 
in this model). The final abundance is 12 + log (O/H) ^ 6.4, ~ 
0.3 dex less than model ICSV. Although in this model also the 
number of stars in the interval 4-7 M (i.e. the main produc- 
ers of primary nitrogen) is reduced, the final log (N/O) is much 
larger (~ 0.4 - 0.5 dex) compared to model ICSV. Finally, due 
to the fact that both massive and low-mass stars contribute to 
the final carbon abundance, the final log (C/O) does not deviate 



much from the value found in model ICSV, being only ~ 0. 1 dex 
larger. 

5.3. Model NCSM 

Model NCSM has an initial setup aimed at reproducing the gross 
characteristics of NGC1569. For this model we assume, accord- 
ing to the work of Angeretti et al. (2005) three episodes of SF: 
a most recent one occurred between 37 and 13 Myr ago, at a 
rate of 0.13 M yr _1 ; an intermediate episode, commencing 150 
Myr ago and finishing 40 Myr ago, at a rate of 0.04 M yr~' and 
an older episode of SF, ending 300 Myr ago (therefore implying 
150 Myr of inactivity between this episode and the intermediate 
one) and commencing 600 Myr ago, at a rate of 0.05 M yr _1 . 
The setup is the same as the model NGC - 5 described in Paper 
II, namely the total mass inside the galaxy is 1.8 x 10 8 M , but 
we add continuously created clouds. At variance with the models 
ICSV and ICKV, the clouds are created every 5 Myr, therefore 
their masses are between 2 and 6.5 x 10 5 M , not far from the 
values actually observed in the complex of H i clouds spiraling 
around NGC 1569 (Miihle et al. [2005) and significantly larger 
than the ones considered in the previous model with continu- 
ous creation of clouds. Moreover, due to the assumption made 
in Paper II about the nucleosynthetic prescriptions, also in this 
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case we consider yields (from both massive and IMS) taken from 
Meynet & Maeder J2002b . 

5.3.1. Dynamical results 

The dynamical evolution of this model is affected by the larger 
masses of the clouds and by their lower creation rate. It is there- 
fore less frequent the formation of clouds close in space and 
time, able therefore to mutually interact. The evolution of the 
model in the first ~ 100 Myr is shown in Fig. QT| The impact 
of these clouds on the development of large-scale outflows is 
stronger than in the above-presented cases. We can clearly see 
the collision of a cloud with the expanding supershell in the 
lower right panel of this figure. These massive clouds lead also 
to a stronger evaporation rate, therefore they affect significantly 
the energy budget of the model. The total thermal energy is ~ 
35 - 40% smaller than the one attained by model NGC - 5 (de- 
pending on the evolutionary time). Again, as we have noticed 
in Sect. 15.1.11 in spite of the reduced superbubble power, the 
supershell-cloud interaction can create holes in the supershell 
from which highly pressurized gas can escape. This is visible 
for instance in the central right and in the upper central panel of 
Fig-Qj] At ~ 100 Myr the dimension of the superbubble is only 
slightly smaller than the one attained at the same evolutionary 
time by model NGC - 5. This demonstrates once again that the 
superbubble luminosity is just one factor, in many cases not the 
leading one, in determining the development and the shape of 
large-scale outflows. Very important factors are also the density 
distribution and the interaction with clouds. 

5.3.2. Chemical results 

Similarly to what is seen in Sect. 15.1.21 the dilution effect of 
the clouds overcomes the "cap" effect, due to the delayed devel- 
opment of an outflow. Therefore, model NCSM shows a lower 
metallicity compared to model NGC - 5. In particular the oxy- 
gen abundance, after an initial phase of ~ 250 Myr in which it 
grows constantly, remains in the range 12 + log(0/H) ~ 7.6 - 
7.7, therefore 0.3 - 0.4 dex below the value attained by model 
NGC - 5. Similarly to the diffuse model, the last intense burst 
of SF has some effect on the chemical evolution of this model, 
resulting in an increase of ~ 0. 1 - 0.2 dex. Also the log(N/0) 
is reduced by ~ 0.2 - 0.3 dex compared to the corresponding 
diffuse model. Since we had noticed in Paper II that the results 
model NGC - 5 matched well the observed chemical composi- 
tion of NGC 1569 (taken from Kobulnicky & Skillman |T997l ), we 
can point out that the inclusion of a cloud complex worsens the 
agreement between model results and observations. In order to 
match the observed chemical composition, one should therefore 
reduce the total initial mass, in order to diminish the gas fraction 
and therefore increase the metallicity. Playing with this param- 
eter is allowed by the present uncertainties about the total mass 
of NGC1569 (Stil & Israel |2002l Miihle et al. 15005V but, as al- 
ready pointed out in the Introduction, our main focus is not the 
quest for the best setup able to reproduce the chemical compo- 
sition of specific objects, but the study of the effect of a cloud 
complex and the differences with a model in which the gaseous 
distribution is smooth. 



6. Discussion and conclusions 

In this paper we have computed the chemical and dynamical evo- 
lution of model galaxies, with structural parameters similar to 



IZwl8 and NGC 1569, but in which a complex of clouds has 
been added, both perturbing the initial gaseous distribution and 
creating clouds, at a rate which equals the SF rate, and with infall 
velocity of 10 km s along the polar direction. The main focus 
of our work has been the comparison of these models with those 
presented in previous publications, in which similar setups but a 
smooth distribution of gas was considered. 

We have seen that the clouds are subject to a variety of dis- 
ruptive phenomena like evaporation (when embedded in a hot 
medium), formation of shocks, development of thermal insta- 
bilities (in particular the Kelvin-Helmholtz instability) and ex- 
pansion due to the larger pressure compared to the surround- 
ing interstellar medium. The average lifetime of the clouds is 
therefore relatively short, depending on the cloud size (which is 
not constant in our simulations) but being of the order of a few 
tens of Myr. In spite of their transient nature, the clouds leave 
a significant imprint on the dynamical and chemical evolution 
of dwarf galaxies. The clouds, when they evaporate inside the 
superbubble, produce mass loading, increase the mean density 
of the cavity and, therefore, enhance the radiative losses (which 
are proportional to the square of the density). This results in a 
significant decrease of the total thermal energy (of the order of 
~ 20 - 40% compared to the diffuse models, depending on the 
assumptions), therefore less energy to drive the development of 
a large-scale outflow. 

On the other hand, the relative motion of supershell and 
clouds, in particular when the clouds infall motion is consid- 
ered, can structure, pierce and create holes and fingers in the 
expanding supershell. These holes destroy the spherical sym- 
metry initially present and favor the rushing out of the highly 
pressurized gas contained in the cavity. Therefore, in spite of the 
reduced thermal energy budget, the creation of large-scale out- 
flows is not suppressed but, in most of the explored cases, only 
slightly delayed. Complex structures and fingers are indeed rela- 
tively common features in galaxies showing large-scale outflows 
like NGC1800 (Hunter fl99%l ) NGC 4214 (MacKenty et al. l2000b 
or NGC1705 (Heckman et al. 2001). The pressure inside the cav- 
ity is reduced compared to diffuse models, therefore in any case 
the total amount of ejected pristine gas is very small (smaller 
than in the models with smooth gas distribution) and, when av- 
eraging the size of the supershell in any direction, it turns out 
to be smaller than in diffuse models. But the piercing of the su- 
pershell can lead to an ejection efficiency of freshly produced 
metals as high as the one attained by diffuse models. 

This has, of course, important consequences on the chemical 
evolution of these objects. Since the differential winds are not 
suppressed, the diminished thermal energy of these models does 
not imply an increase of metals inside the galactic regions. On 
the other hand, the dilution effect of clouds plays a dominant role 
in determining the final metallicity of our model galaxies. Since 
the clouds have primordial chemical composition, their destruc- 
tion and mixing with the surrounding medium reduces the to- 
tal chemical composition without altering the abundance ratios. 
This produces a final metallicity ~ 0.2 - 0.4 dex smaller than the 
corresponding diffuse models. 

We have examined the effect of a different choice of the 
IMF slope and of the nucleosynthetic set of yields (in partic- 
ular for what concerns intermediate-mass stars). Flatter-than- 
Salpeter IMF slopes lead to an excessive production of energy, 
able to unbind most of the gas before the end of the simula- 
tion. On the other hand, in models with steeper IMF the devel- 
opment of large-scale outflows is almost completely suppressed. 
Different sets of intermediate-mass stars yields affect in partic- 
ular the log(N/0) ratio. Renzini & Voli ( 119811 1 yields tend to 
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Fig. 11. Same as Fig. [I0]but for model NCSM. For reference, the density contours of model NGC - 5 (presented in Paper II) are 
shown in the upper right panel. 



overestimate the primary production of nitrogen. When com- 
pared to the results of models implementing van den Hoek & 
Groenewegen (1997 ) yields, the results differ by ~ 0.3 dex. Due 
to the assumption of a metallicity-dependent cooling function, 
also the dynamics is affected by the choice of the nucleosyn- 
thetic prescriptions. 

Our main results can be briefly summarized as follows: 

- the clouds suffer thermal instabilities, formation of shocks 
and evaporation, therefore their lifetimes is limited to a few 
tens of Myr. 

- In spite of that, they are able to increase the main density of 
the cavity, provoking a reduction of the total thermal energy 
by ~ 20 - 40% compared with a diffuse model. 

- The interaction clouds-supershell leads to strong structuring 
and piercing of the shell (in particular for models with con- 
tinuous creation of infalling clouds), allowing the venting out 
of metals in spite of the reduced thermal energy. The devel- 
opment of large-scale outflows is therefore generally delayed 
but the ejection efficiency of metals remains unchanged. 

- From a chemical point of view, the effect of the clouds is 
to significantly reduce the total metallicity of the galaxies, 
without altering the abundance ratios. 
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